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Abstract 

The theoretical treatment of alpha decay by Gamow is revisited by investigating the quantum 
time scales in tunneling. The time spent by an alpha particle in front of the barrier and traversing 
it before escape is evaluated using microscopic alpha nucleus potentials. The half-life of a nucleus is 
shown to correspond to the time spent by the alpha knocking in front of the barrier. Calculations for 
medium and super heavy nuclei show that from a multitude of available tunneling time definitions, 
the transmission dwell time gives the bulk of the lifetime of the decaying state, in most cases. 

PACS numbers: 03.65.Xp, 21.10.Tg, 23.60.+e 



1 



The concept of quantum tunneling has been with us for decades and has been successfully 
applied in many branches of physics. It essentially expresses the fact that in quantum 
mechanics there is always a non-zero finite probability for a particle to go from one region 
to another even if the regions are separated by a large but finite potential barrier and the 
particle carries a kinetic energy less than the height of the barrier. Some of the very first 
applications were by Gamow [l| and by Gurney and Condon to study the alpha decay 
of radioactive nuclei. It was however soon noticed that the tunneling phenomenon was 
not restricted to nuclear physics, but was rather a general result of quantum mechanics 
which is now commonly used to study the physics of semiconductors and superconductors, 
in constructing electron tunneling microscopes, to find the lifetimes of the newly discovered 
super heavy elements and sometimes even to understand the early cosmology of the 
universe Q|. Its many ramifications reached the realm of atomic physics [5] with ultracold 
atoms [6| , as well as chemistry [7|, l8| and biology . In spite of the fact that tunneling 
is now widely applied in many fields, the question of how much time does the particle 
require to tunnel and what connection does it have with the physically measured lifetimes 
of spontaneously decaying objects still remains unanswered in the cases where we encounter 
a bound state before tunneling. 

In the present work we attempt to answer the above question. To be specific, within the 
framework of semi-classical approximations, we derive expressions for the times spent by a 
tunneling particle in front of the barrier as well as within the barrier. Examining the known 
expressions for the lifetimes of decaying (metastable) states (or resonances), a connection is 
found between the so-called dwell times in tunneling and the lifetimes of metastable states. 
The dwell time is considered to be a measure of the average time spent by a particle in 
a given region of space. The concept was first introduced by Smith 10|] in the context of 
quantum collisions and to derive a lifetime matrix for multichannel resonances. In the one- 

n 

dimensional case, it was first introduced by Biittiker Indeed, earlier in 1966, Baz had 
proposed 12|] the use of Larmor precession as a clock to measure the duration of quantum 
mechanical collision events and Rybachenko 13[ had applied the method to the simpler case 
of particles in one dimension. The work of Biittiker was an extension of these works. Biittiker 
and Landauer also defined a traversal or interaction time of transmitted particles in 
tunneling, which, as we shall see later also finds a new physical significance in the tunneling 
process. More recently, the dwell time formalism for the transition from a quasilevel to a 
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continuum of states was discussed in the context of electron and alpha particle tunneling by 



P. J. Price 



15| . Some other recent applications of dwell time can be found in [l6l |. 



In what follows we shall present the formalism connecting the dwell times with the decay 
widths (and hence the half-lives) and conclude with the help of realistic examples that the 
dwell times of transmitted alpha particles in the region before entering the barrier correspond 
to the half-lives of radioactive nuclei decaying by alpha decay. Apart from the dwell time, 
there exist other tunneling time concepts in literature and the subject in principle is replete 
with controversies 17|]. The controversy stems from definitions of phase times or the so- 



called "group delay times" which can also predict superluminal tunneling velocities. The 
Hartman effect jl?!] for example, states that the tunneling time becomes independent of 
the thickness of the barrier length for thick enough barriers and thus results in unbounded 
tunneling velocities. Such controversies, however, are not relevant for the calculations done 
in the present work. We do not evaluate the phase time but rather the dwell time and the 
barrier region itself makes a small contribution to the total dwell time (with the bulk coming 
from the region in front of the barrier). Hence we do not expect saturation effects such as 
the Hartman effect resulting in interpretations of superluminal alpha propagation. 

Though the nuclear potentials used in this work are not the most modern ones, these 
potentials and the semiclassical approach [18| used suffice to illustrate the main findings of 
this work. There exist potentials for which the semiclassical WKB approach may not be 
appropriate. However, it is known to work reasonably well for the alpha tunneling problem. 
This approach is indeed commonly used in literature for the evaluation of the standard 
Gamow factor and prediction of half-lives of heavy and super heavy nuclei [lol . We shall see 
below that the dwell time (which one would generally not expect to be an experimentally 
measurable quantity) in the region in front of the barrier corresponds to the commonly used 
definition of the measured half-life. Since the WKB approximation is widely used for the 
calculation of half-lives, we too use the WKB wave function for the evaluation of dwell times 
in the present work. 

For an arbitrary barrier V{x) in one-dimension (a framework which is also suitable for 
spherically symmetric problems), confined to an interval {xi,X2), the dwell time is given by 
the number of particles in the region divided by the incident fiux j: 
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Here is the time independent solution of the Schrodinger equation in the given region. 
Though we shall restrict to using a semiclassical approximation for the wave function in the 
present work, one can always define a dwell time as above in a given region of space, be it 
with an exact or an approximate wave function. The standard definition of dwell time is 
the time spent in the region (xi , X2) regardless of how the particle escaped (by reflection 



or transmission) and j = hko /jj, [20] (with ko = y/2jlE/fi) for a free particle. However, 
one can also define transmission and reflection dwell times for the particular cases when the 
particle is bound in a region and later either got transmitted or reflected. The flux j in 
these cases would get replaced by the transmitted or reflected fluxes, jr = hko\T\'^/fi and 
j/j = hkolRl'^/n I2II respectively. The current Jt is the particle's flux in the region HI (see 



Fig. 1). One would then obtain 211] 
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where |Tp and |-Rp are the transmission and reflection coefficients (with |Tp + = 1 due 
to conservation of probability) and T£, t = J \ '^\'^dx/jT and T£, ji = J I'^l'^dx/jn, define the 
transmission and reflection dwell times respectively. The traversal time defined by Biittiker 
is somewhat different and is given as, 

= /r nW) 

where, k{x) = ^2fi{\V{x) — E\)/h with E being the kinetic energy of the tunneling particle 
and /i the reduced mass. Having defined the tunneling times relevant to the present work, 
we shall now apply them to the study of the alpha decay of nuclei. 

We study the alpha decay of nuclei as a tunneling of the a through the potential barrier of 
the alpha-daughter nucleus system using a semi-classical approach. Typically, one considers 
the tunneling of the a through a spherically symmetric r-space potential (see Fig. 1) of the 
form, V{r) = Vn{r) + Vc{r) + {I + 1/2)^/2/2 r^, where Ki(r) and Vc(r) are the attractive 
nuclear and repulsive Coulomb parts of the a-(daughter) nucleus potential, r the distance 
between the centres of mass of the daughter nucleus and alpha and /x their reduced mass. 
The last term represents the Langer modified centrifugal barrier 22|. Writing the wave 



functions in region I and region II (up to a normalization factor) using the semi-classical 



Wentzel-Kramers-Brillouin (WKB) approximation 23] 

'k{r 
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FIG. 1: Typical potential y(r) in an alpha-nucleus tunneling problem, ri, r2 and rs are the 
classical turning points for a given kinetic energy E of the tunneling particle. 
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where k{r) = ^2fi{E — V{r))/h and ^(r) = ^2fi {V{r) — E)/h, the dwell times and 
in regions I and II respectively are given as follows: 
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where the second line follows from replacing the squared cosine term by (1/2). The dwell 
time in region II is, 

rrs dr 



exp 



-2 



dr' n{r') 



(6) 



h Jr2 K{r) 

Though the semiclassical WKB approach which is still often used for the evaluation of decay 



widths 



24( 1 is sufficient for the objectives of the present work, better methods such as the 



Gamow-state formalism [25|] do exist. In fact, earlier, starting in the seventies, the general 
formula for the lifetime of a nucleus decaying by a-decay was obtained on the basis of a 
Gamow-state formalism in (26|]. If we now consider the standard definition of the WKB 

h"^ r dr 



decay width [27|, 



r{E) = p. 



2/i 



k{r) 



-2 K{r)dr 
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{Pa is the pre-formation probability of the alpha cluster inside the radioactive nucleus de- 
caying later by alpha particle emission) and compare it with equation ([5]), taken along with 
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the fact that the transmission coefficient |T| 
result of this work, 



-2 K(r)dr 



we arrive at the first main 



T{E) = P^n\T\' [rh{E)]-' = P^n [T^^TiET'- 



(8) 



The decay width is given by the inverse of the transmission dwell time in the region in front 
of the barrier. This imphes that the half-hfe (which is evaluated at the energy E = Q, where 
Q is the amount of energy released in the decay), 



^ln2 



In 2 



n/2 = = ^D,T) (9) 

is essentially given by the transmission dwell time in region I. It is of further interest to note 
that the frequency of assaults at the barrier, can be written as the inverse of the time 
required to traverse the distance back and forth between the turning points ri and r2 as 



28|, 
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However, from equation ([5]), it follows that v = [r^] ^ and taken along with equation ([8]) 
with Pq, = 1, we can see that the number of assaults that the particle makes before tunneling 



is Na 



D,T 



m 



In Table I, we list the transmission dwell times for four medium heavy nuclei with / = 
and two recently studied super heavy nuclei In the absence of much information on 

the / values of the super heavy nuclei studied, we assume the angular momentum / = 0. 
The case of the light nucleus ^Be (with / = 2) is presented separately in Table II. The 



alpha-alpha potential for ®Be is described analytically and is taken from 



30|. The alpha- 



nucleus potential for the medium and super heavy nuclei is constructed using a double 
folding model with realistic nuc 



used in some recent works 
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.eon-nucleon (A^iV)interactions as given in |30l] and also 
3l| . An accurate estimate would however demand the latest 
available A^A^ potentials [3^ and considerations of the four-particle correlations in nuclei 33|. 
The Coulomb potential is also obtained via a double folding procedure where the matter 
densities of nuclei are replaced by their charge densities. The details of the potentials used 
in the present work can be found in 34 1. 

We also list the transmission dwell times in region II, in the tables. To evaluate the total 
time spent in regions I and II, one would start with the definition of the dwell time as, 

m \\^(r)\'^dr 
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Jri 
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(11) 



The lifetime of the decaying nucleus should in principle be related to the total transmission 
dwell time spent by the a in the two regions and not just in region I as given by ([9]). Defining 
such a time as, say, 

r^^'(^)=^L,T(^) + rD,T(^), (12) 

the half-life expression in ([9]) would rather take the form. 

In 2 



^fuU 



'1/2 - pfuU '^D,T^ 

a 

with P^"'' being the preformation factor in this case. For the medium heavy nuclei, agreement 
with experimental half-lives is obtained with P^"^^ ~ 0.4 — 0.55. These values are close to 
those obtained in literature [l9| with similar potentials. The P^^" values for the super heavy 
ones are around 0.2, which could change with possible higher angular momenta for these 
nuclei. In principle, the values of P^"^^ depend on the details of the nuclear potentials which 
in turn affect the tunneling probabilities. The values of P^" ^' found here should hence not be 
interpreted as a constraint on four-particle correlations 33|| in microscopic nuclear structure 
models. A detailed study of the preformation factors for super heavy alpha emitters can be 
found in [35|. The half- lives are evaluated at the experimental Q values, i.e., sX E = Q m. 
equation 



full 



(13) 



TABLE I: Comparison of the calculated half lives (from the transmission dwell times) with exper- 
iment, for medium and super heavy alpha emitters. 



Radioactive 
Nucleus 


Q value 
(MeV) 


ln2 t\) rp 

(s) 


ln2 rj^j. 
(s) 


ln2 {t^^t + ^d,t) 
(s) 


ri/2 (expt) 
(s) 


pfuU a 


Number of 
assaults 


lO^Te 


3.445 


1.506 


0.261 


1.767 


4.286 


0.41 


7 X 10^1 


169t^ 


6.151 


0.466 


0.069 


0.535 


1.28 


0.42 


2 X 10^1 


173 An 
79 


6.836 


9.29 X 10-3 


1.38 X 10-3 


10.67 X IQ-^ 


26.6 X 10-3 


0.40 


4 X 10^9 




2.508 


2.69 X 10^5 


3.74 X 10^4 


3.06 X 10^5 


5.68 X 10^5 


0.54 


1 X 10^7 


273t-)„ 


11.368 


2.038 X 10"^ 


2.76 X IQ-^ 


2.32 X 10"^ 


17 X 10"^ 


0.14 


9 X 10^^ 


277112 


11.3 


1.16 X IQ-^ 


1.52 X 10-^ 


1.31 X 10"^ 


6.90 X 10-^ 


0.19 


5 X lO^'^ 



ln2 (tI,^^ + tII_^) / Ti/2(expt) = 



full 



We find that the major bulk of the half-life of a medium or super heavy radioactive 
nucleus is spent in region I, in front of the barrier before tunneling. Though the time spent 
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TABLE II: Same as Table I but for the light nucleus ®Be. 



Radioactive 
Nucleus 


Q value 
(MeV) 


ln2 Tjj rp 

(s) 


(s) 


ln2 {t}j^t + t^^t) 
(s) 


n/2 (expt) 
(s) 


pfull a 

a 


Number of 
assaults 


^Be (1=2) 


3.1218 


3.7 X 10-22 


2.56 X 10-22 


6.26 X 10-22 


3.02 X 10-22 


2.1 


1 






















ln2T|) 


ln2 


ln2 {rj, + rlj) 












1.59 X 10-22 


2.19 X 10-22 


3.78 X 10-22 


3.02 X 10-22 


1.25 


1 



ln2 {t},^^ + t1,V) / Ti/2(expt) = pf' 



full 



inside the barrier is much smaller, it is not negligible and should be added to that in region 
I to obtain the total time spent in tunneling. The number of knocks (assaults) made at the 
barrier is inversely proportional to the transmission coefficient which for the heavy nuclei is 
extremely small leading to a huge number of assaults by the alpha at the barrier. With |i?|2 
being almost unity in these cases, the reflection dwell times, td^h as well as the average dwell 
times td are orders of magnitude smaller as compared to td^j- as well as the half-lives. The 
times spent in region II (the barrier) are always an order of magnitude smaller. The case of 
^Be (with / = 2 studied here) is however, very different. With the transmission coefficient 
being reasonably big (|T|2 = 0.86 at the Q value of 3.1218 MeV), the times, r^ji, t^ t and 
td are all comparable and of the same order of magnitude in both regions. It makes more 
sense then, to compare the standard average dwell time definition with the experimental 
half-life. An interesting implication of the large transmission coefficient here is that the 
alpha in ^Be escapes after one knock. 

In conclusion, we can say that the present work provides a new look at the physics of 
the quantum tunneling problem in general and the alpha decay problem in particular. In 
principle, one does not really know how alpha particles or clusters of heavier nuclei are 
formed inside the nucleus. In fact, there are arguments that they can be formed in virtual 



states [36|] different from the states of those objects in free space. Such objects could undergo 
some restructuring in the nucleus [s^] which can then not be described by a simple potential 
consideration. The present work does not attempt to address these issues. The conclusions 
drawn in this work are hence limited to the picture of a preformed cluster (say the alpha) 
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inside the nucleus, which tries to tunnel through a potential barrier. Further, we note that 
in contrast to the lifetime definition for the decay of an elementary particle, say muon decay, 
the definition of lifetime in the context of a tunneling problem is somewhat different. Here 
the tunneling object has to traverse a region of space, thus giving rise to the question of how 
long does it require for a particle to tunnel (for a general discussion, see 38|]). The dwell 
time concept answers this question. In the alpha tunneling problem studied in the present 
work, it is found that the dwell time of the a particle in front of the barrier as well as within 
the barrier is important. This can be relevant not only to the alpha radioactive nuclei of 
medium mass and the recent discoveries of super heavy nuclei, but also to other branches of 
science and physics where the quantum concepts of time in conjunction with bound states 
and tunneling might prove fruitful to have a measure of the speed of a reaction. 
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